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Moments (Variance (a^), Skewness(5'), Kurtosis(K)) of multiplicity distributions of conserved 
quantities, such as net-baryon,net-cliarge and net-strangeness, are predicted to be sensitive to the 
correlation length of the system and connected to the thermodynamic susceptibilities computed 
in Lattice QCD and Hadron Resonance Gas (HRG) model. In this paper, we point out several 
measurement artifacts like centrality bin width effect, auto-correlation and centrality resolution effect 
in the moment analysis of net-proton multiplicity distributions in high energy heavy-ion collisions 
using the UrQMD model. We present methods to overcome these artifacts so that the extracted 
moments could be used to obtain physical conclusions. In addition we present methods to properly 
estimate the statistical errors in moment analysis. 

PACS numbers: 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 

One of the main goals to perform heavy-ion colli- 
sion experiments is to explore the phase structure of 
the QCD matter, such as search for the QCD critical 
point [l|, Q and the phase transition boundary between 
hadronic matter and partonic matter. Finite tempera- 
ture Lattice QCD calculations predict that the transi- 
tion from nuclear matter phase to Quark Gluon Plasma 
(QGP) phase is a smooth cross-over at high temperature 
and small ub 01i ^^'^ ^ ^^^^ order phase transition at 
large fj,B (j]. The corresponding transition temperature 
is about 170~190 MeV Q. Although many efforts have 
been made by theorist and experimentalist to locate the 
CP (Critical Point), its location or even existence is still 
not confirmed yet [g. The first principle Lattice QCD 
calculation at finite fiB is difhcult to carry out due to the 
sign problem P, Q. Although many techniques, such as 
re- weighting, image baryon chemical potential and Tay- 
lor expansion, have been developed to overcome those 
problems and make the numerical QCD calculations at 
finite fiB region, but with large uncertainties 0. 

In heavy- ion collision, the moments (Variance (cr^), 
Skewness(S'), Kurtosis(K)) of distributions of conserved 
quantities, such as net-baryon, net-charge and net- 
strangenesSjare predicted to be sensitive to the correla- 
tion length 3 and to be connected to the thermodyiiamic 
susceptibilities computed in Lattice QCD [l|, |^, [l^j and 
in the Hadron Resonance Gas (HRG) [ll| model. QCD 
based model calculation show that moments of particle 



multiplicity distributions are proportional to the higher 
power of the correlation length, such as third order cu- 
mulant < {6N)^ >^ ^^-^ and fourth order cumulant < 
{6Ny > -3((5iV)2 In statistics, skewness and 

kurtosis are widely used to characterize the properties of 
probability distributions. Skewness is used to describe 
the asymmetry property of distributions, while kurtosis 
describes their pcakness. The higher moments of baryon 
number distributions are also directly connected to the 
corresponding baryon number susceptibilities in Lattice 
QCD and HRG model calculations [HI, for e.g. the 

/-ON 

third order susceptibility of baryon number {Xb ) 
lated to the third moment (< (SNb)^ >) of baryon num- 
ber distributions as Xb^ =< {SNbT > /VT^; and V,T 
are volume and temperature of system respectively. As 
volume of the system is hard to determine, the suscep- 



tibility ratio, such as Xb /Xb ^^'^ Xb' /Xb'' used 
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and xL^Vxg\ 
to compare with the experimental results. It has been 
found that the baryon number susceptibility ratio can 
be directly connected to the moment products of net- 
baryon number distributions as Ka^ = X^Vxb'' and 

^'^ = XbVxb^ 0- Theoretical calculations also demon- 
strate that the experimental measurable net-proton num- 
ber Np-p = Np — Np (number of protons minus number 
of anti-protons) event-by-event fluctuations can reflect 
the baryon, charge and strangeness number fluctuations 
and also sensitive to the correlation length [l^. Thus, 
we can use the higher order moments of the net-proton 
multiplicity distributions to search for the QCD critical 
point 

The paper is organized as follows. In section II we 
briefly discuss the UrQMD model used for the analysis 
presented here. In section III wc define the observables. 



2 



Sections IV, V and VI discuss the measurement artifacts, 
centrahty bin width effect (CBWE), centrality resolution 
effect (CRE) and Auto-correlations effect (ACE), respec- 
tively. In section VII we present different methods of esti- 
mating the statistical errors for moment analysis. Finally 
in section VIII we present a summary of the work. 



II. URQMD MODEL 

The Ultra Rclativistic Quantum Molecular Dynamics 
(UrQMD) is a microscopic many-body approach to p+p, 
p+A, and A -I- A interactions at relativistic energies and is 
based on the covariant propagation of color strings, con- 
stituent quarks, and diquarks accompanied by mesonic 
and baryonic degrees of freedom. Furthermore it in- 
cludes rescattering of particles, the excitation and frag- 
mentation of color strings, and the formation and decay 
of hadronic resonances. UrQMD is a transport model 
for simulating heavy-ion collisions in the energy range 
from SIS to RHIC (even in LHC). It combines different 
reaction mechanism, and can provide theory simulation 
results of various experimental observables. The main 
parts of UrQMD model are: two body reaction cross 
section, two body potential and decay width. More de- 
tails about the UrQMD model can be found in the refer- 
ence [ist . 

We performed our calculations with UrQMD model 
for Au-I-Au collisions at ^/s^= 7.7, 11.5, 19.6, 27, 39, 
62.4, 200 GeV and the corresponding event statistics 
are 107, 108, 50, 65, 93, 29 and 28 million, respec- 
tively. We measure event-by-event wise Np_p at the mid- 
rapidity (|y| < 0.5) and within the transverse momentum 
0.4 < pt < 0.8 GeV/c. This is same kinematic range 
used in moment analysis of net-proton distributions as 
in the STAR experiment [H [l6|-[lil . 



III. OBSERVABLES 

In statistics, probability density distribution functions 
can be characterized by the various moments, such as 
mean (M), variance (cr^), skcwncss (5) and kurtosis (k). 
Before introducing the above moments used in our anal- 
ysis, we would like to define cumulants, which are alter- 
native methods to the moments of a distribution. The 
cumulants determine the moments in the sense that any 
two probability distributions whose cumulants are iden- 
tical will have identical moments as well, and similarly 
the moments determine the cumulants. 

In the following, we use N to represent the net-proton 
number Np-p in one event. The average value over whole 
event ensemble is denoted by < >. where the single 
angle brackets are used to indicate ensemble average of 
a event-by-event distribution. 

The deviation of N from its mean value is defined by 



The various order cumulants of evcnt-by-evcnt distri- 
butions of a variable N are defined as: 



C2,N 

Can 



<N> (2) 

< i6N'f > (3) 

< {SNf > (4) 

< {5N)^ > -3 < {SN)^ >2 (5) 



An important property of the cumulants is their ad- 
ditivity for independent variables. If X and Y are two 
independent random variables, then we have Ci^x+Y = 
Ci^x + Ci.Y for zth order cumulant. This property will 
be used in our study. 

Once we have the definition of cumulants, various mo- 
ments can be obtained as: 



M = Ci^AT, cr^ = C2,N,S — 



a 



3,N 



a 



13/2' 



K = 



4,Af 



iC2,Ny 



(6) 



And also, the moments product kct^ and Sa can be 
expressed in terms of the ratio of cumulants like: 
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, Sa 
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C2. 
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(7) 



With above definition of various moments, we can cal- 
culate various moments and moment products from the 
measured event-by-cvent net-proton distribution for each 
centrality. 



IV. CENTRALITY BIN WIDTH EFFECT 
(CBWE) AND CORRECTION 

To avoid the ACE and CRE (will be discuss later), 
the centralities in the UrQMD calculations are defined 
by the charged kaon and pion multiplicity within \r]\ < 2. 
Experimentally, the centrality is usually expressed as a 
percentage of the total cross-section, such as 0—5% (most 
central), 30 — 50% (semi- peripheral), which indicates the 
fraction of a data sample (corrected for efficiency) relative 
to all possible collision geometries (impact parameters). 
Before calculating various moments of net-proton distri- 
butions for one wide centrality bin, such as 0—5%, 5—10% 
, we should consider the so called Centrality Bin Width 
Effect (CBWE) arising from the impact parameter (or 
volume) variations due to the finite centrality bin. This 
effect must be eliminated or at least reduced, as an ar- 
tificial centrality dependence could be introduced due to 
finite centrality bin width. 

To formulate and demonstrate the centrality bin width 
effect, we write the event-by-event net-proton distribu- 
tions in one centrality as 



SN ^ N- < N > 



(1) 



(8) 
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FIG. 1: (Color online) The centrality dependence of the mo- 
ments products (Sa) of net-proton multiplicity distributions 
for Au-FAu coUisions at ^s^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. The dots represent the result of 
32 centralities. 
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FIG. 2: (Color online) The centrality dependence of the mo- 
ments products (ko-^) of net-proton multiplicity distributions 
for Au+Au collisions at ^s;7;7=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. The dots represent the result of 
32 centralities. 



, where the and f^^\N) are the weighted and net- 
proton distributions for ith impact parameter in one cen- 
trality, respectively. From eqs. (2)- (5) and (8), we can 
calculate the various order cumulants for the distribution 
P{N) as below: 



C2,N 
C3,N 

Can 



= ^ uj.C'Ij^ ^ < TV >i (9) 

i i 

- (E^'^2,w) + ^2,c;,„ (10) 

i 
i 
i 



3 X (C2 (^j + 3 ^ 



2,Cl 



(12) 



, where the — li2,3,4) are the fc order 

cumulant for net-proton distribution f^^\N)] the C^, 
[X ~ C'mTV' '^'"^ = 1;2,3,4) are k*-^ order cutmulant 



X 



for random variable X = 
abihty distribution Prob{X) 
>< Y > {X 



^ under the discrete prob- 

=Wi; the X I y~ 
- < X >< Y > {X = Clj^, k = 1,2,3,4) are the 
first order joint cumulant for random variable X, Y un- 
der the discrete probability distribution Prob{X, Y)=uji. 
We may find that the higher order cumulants Ck,N 
{k = 1,2,3,4) can be expressed by the addition of two 
parts, one is the weighted average of the same order cu- 
mulant of each sub-distribution /'^'^(A^), and the other 
part is the cumulant of lower order under the discrete 
weighted distributions, which stems from the fluctuation 
of impact parameters within the centrality and results in 
the centrality bin width effect. 

Experimentally, the centrality is usually determined by 
an observable like the charged particle multiplicity and 
the smallest centrality bin is determined by a single mul- 
tiplicity value. Generally, experiments report the results 
for a wider centrality bin, such as — 5%, 5 — 10%, ...etc., 
to supress the statistical fluctuations. To eliminate the 
centrality bin width effect, we have developed a tech- 
nique, called Centrality Bin Width Correction (CBWC), 
to calculate the various moment for each multiplicity in 
one wide centrality bin. The method has the weighted 
averaged by the number of events in each multiplicity. 



r 

r 

r 

r 

r 

J2nr 



r 
r 



(13) 



(14) 



(15) 
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, where the is the number of events in r*'* multiphc- 
ity and the corresponding weight ujr = rirf^nr. 

r 

Fig. [T] and Fig. [5] show the centrahty dependence of 
the moment products {Sa,Ka'^) of net-proton multiphc- 
ity distributions for Au+Au colhsions ^s„j^=7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV in the UrQMD model. 
The open circle and open cross in Fig. [1] and [2] repre- 
sent the results calculated with and without technique 
of CBWC in the nine centralities (0-5%, 5-10%, 10-20%, 
20-30% ...), respectively. The solid circles show the result 
of centrality dependence of 32 bins in centrality (0-2.5%, 
2.5-5%, 5-7.5%...), in which the CBWC is not applied. 
We can see clearly that the results with CBWC are very 
different from the results without CBWC. It means the 
CBWE do have a significant contribution to the value of 
the moment analysis and the CBWC will make the value 
of the moment products smaller by reducing the eS'ect of 
variation of volume in one wide centrality bin. On the 
other hand, the results of 32 centralities, for which the 
CBWE is small, are good agreement with the result from 
nine centralities with CBWC, which indicates the effec- 
tiveness of the proposed CBWC method. So, before we 
perform higher order moments calculation in one wide 
centrality in heavy-ion collisions, the CBWC should be 
applied to reduce the effect of variation of the volume 
within the centrality bin width. 

We also try to use error weighted method to do the 
CBWC (open star in Fig. [1] and [5]) by replacing the weight 
factor rir in Equ. (13)-(15) with 1/error^ for each single 
multiplicity. The error is obtained with Delta theorem, 
which will be discussed later. We find that the results 
with error weighted CBWC are consistent with the re- 
sults with events weighted CBWC and the results with 
32 centralities for Sa, but not for kct^. It means the error 
weighted method can not be used for CBWC, which may 
be due to the error being not only related to the number 
of events but also the moments itself. 
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FIG. 3: (Color online) The centrality dependence of the cr^/M 
of net-proton multiplicity distributions for Au-|-Au collisions 
at ^/s]^=7.7, 11.5, 19.6, 27, 39, 62.4, 200GeV in UrQMD 
model. 
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V. CENTRALITY RESOLUTION EFFECT 

The definition of the centralities in heavy-ion collisions 
is not unique and can be defined through different quan- 
tities. A frequently used quantity is the so called im- 
pact parameter b, defined as the distance between the 
geometrical centers of the colliding nuclei in the plane 
transverse to their direction. Other variables, such as the 
number of participant nucleons, Npart and the number of 
binary collisions, Ncoii, can be also used. Unfortunately, 
those geometry variables can't be directly measured in 
the heavy-ion collision experiment and the collision cen- 
trality is determined from a comparison between exper- 
imental measures such as the particle multiplicity and 
Glauber Monte-Carlo simulations [l^. Particle multi- 
plicity, not only depends on the physics processes, but 
also reflects the initial geometry of the heavy-ion colli- 
sion. This indicates that relation between measured par- 



FIG. 4: (Color online) The centrality dependence of the mo- 
ments products Sa of net-proton multiplicity distributions 
for Au-f An collisions at ,/s;;^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. 



tide multiplicity and impact parameter is not one-to-one 
correspond and there are fluctuations in the particle mul- 
tiplicity even for a fixed impact parameter. Thus, it could 
have different initial collision geometry resolution (cen- 
trality resolution) for different centrality definitions with 
particle multiplicity. Fig. [3] shows the centrality depen- 
dence of cr^/Af of number of participant nucleons (Npart) 
distributions for Au-f- An collisions at GeV in 

UrQMD model. As the Npart can reflect the initial geom- 
etry (volume) of the colliding nuclei, the a^/M of Npart 
multiplicity distributions can be regarded as the central- 
ity resolution for a certain centrality definition. It means 
if less particles are used to define the collision centrality. 
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FIG. 5: (Color online) The centrality dependence of the mo- 
ments products Ko^ of net-proton multiplicity distributions 
for Au-HAu collisions at ys^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. 



pion used to determine the centrality. When we increase 
the T] range {\r]\ < 1,1.5,2), the values of Sa and kct^ 
will decreases. We observe significant difference for mo- 
ment products (S'cr, Kcr^) for the different rj range of the 
centrality definition. The behavior can be understood 
as due to reduced contribution from volume fluctuations 
arising from centrality definition. On the other hand, the 
moment products {S<7,Ka'^) show saturation at |?7| < 2. 
The centrality resolution effect by now gets largely sup- 
pressed. Fig. [HI shows the energy dependence of moment 
product {Sa, na^) of net-proton multiplicity distributions 
for three difTerent centralities (0-5%, 30-40%, 70-80%) 
with four different rj coverage in Au-|-Au collisions. We 
can find that the fourth order fluctuation is more sensi- 
tive to the centrality resolution effect than the third order 
fluctuation, and it has larger effect in the peripheral col- 
lision as well as at low energies. Thus, we should use a 
lager 77 coverage in the centrality definition for the real 
experimental moment analysis to reduce the centrality 
resolution effects. 
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FIG. 6: (Color online) The energy dependence of the mo- 
ments products (Sa, Ka^) of net-proton multiplicity distribu- 
tions for Au-I-Au collisions at ^s„„=7.7, 11.5, 19.6, 27, 39, 
62.4, 200GeV in UrQMD model. 
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FIG. 7: (Color online) The centrality dependence of the mo- 
ments products Sa of net-proton multiplicity distributions 
for Au-I-Au collisions at ^^3^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. The open square and circle are 
the results with CBWC and without CBWC at \ri\ < 0.5 for 
centrality definition. The solid square and circle are the re- 
sults with CBWC and without CBWC at I?)! < 2 for centrality 
definition. 



then larger initial geometry (volume) fluctuations will be 
observed. This may affect the moment analysis of parti- 
cle multiplicity distributions. To verify this effect in the 
moment analysis, we use the charged kaon and pion mul- 
tiphcity produced in the final state within \ri\ < 0.5, 1, 1.5 
and 2 to define the centrality in the UrQMD calculations. 
Fig. m and |5| shows the centrality dependence of the mo- 
ment products (Sa, Ka^) of net-proton multiplicity dis- 
tributions for four different 77 range of charged kaon and 



Fig. |7| and ID show the centrality dependence of mo- 
ment products {Scr, na"^) of net-proton multiplicity distri- 
butions for Au+Au collisions with CBWC and without 
CBWC at different centrality resolution {\r]\ < 0.5,2). 
We can see that the CBWE and CRE have larger ef- 
fect in the moment products at low energies than at high 
energies. The magnitude of CBWE in the moment prod- 
ucts is comparable with the magnitude of CRE at low 
energies, while the CBWE has larger effects at higher 
energies. 
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FIG. 8: (Color online) The centrality dependence of the mo- 
ment products Kcr^ of net-proton multiplicity distributions 
for Au-hAu collisions at ^s^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. The open square and circle are 
the results with CBWC and without CBWC at \t]\ < 0.5 for 
centrality definition. The solid square and circle are the re- 
sults with CBWC and without CBWC at |77| < 2 for centrality 
definition. 



auto-correlation effects. 




FIG. 9: (Color online) The proton distribution in the phase 
space of transverse momentum (pr GeV/c) and rapidity 
(y) for Au-|-Au collisions at ^s„„ = 19.6 GeV using the 
UrQMD model. The relation between pr and y is pr = 
mo / \/ sinh^r;/sinh'^y — 1, where the mo is the particle rest 
mass. 



Let's summarize the CBWE and the CRE in the mo- 
ment analysis in heavy-ion collisions. The CBWE stems 
from the variation of the volume when we perform the 
moment calculations within one wide centrality bin and 
the CRE is originated from the initial geometry (volume) 
fluctuations of the colliding nuclei. Both of these two 
effects are resulted from the volume effects and it can 
be effectively suppressed by using the CBWC and more 
particles (larger rj coverage) in centrality definition, re- 
spectively. 



VI. AUTO-CORRELATION EFFECT 

In the above sections, we use the multiplicity of 
charged kaon and pion to define the centrality, not all of 
the charged particles. The reason is to avoid the effect of 
auto-correlation between protons/anti-protons involved 
in our moments analysis and in the centrality definition. 
If the charged particles for the centrality definition con- 
tains the proton and anti-proton used in our moments 
analysis, it will introduce auto-correlation effect. Fig. [5] 
shows the proton distribution in the phase space of trans- 
verse momentum (pr) and rapidity (y) for Au+Au col- 
lisions at ys;;;;=19.6 CcV in the UrQMD model. The 
protons used in the moment analysis are distributed in 
the oblong box (0.4 < pr < 0.8 GeV/c, \y\ < 0.5). The 
solid curve represents {px, y) points with rj = 2. If all 
the charged particles within the rj = 2 curve are used 
to define the centrality, it will definitely introduce the 
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FIG. 10: (Color online) The centrality dependence of the mo- 
ments products Sa of net-proton multiplicity distributions 
for Au+Au collisions at yJT^^T.l, 11.5, 19.6, 27, 39, 62.4, 
200GcV in UrQMD model. 

Fig. [10] and [11] show the centrality dependence of the 
moment products Sa and na^ of net-proton distribu- 
tions at centrality definition with charged particles (open 
squares) and charged kaon and pion (solid circles). An 
obvious difference between the results from the two types 
of centrality definition can be observed. This is because 
that if the protons and anti-protons used in the moments 
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FIG. 11; (Color online) The centrality dependence of the mo- 
ments products Ka^ of net-proton multiplicity distributions 
for Au-HAu collisions at ys^=7.7, 11.5, 19.6, 27, 39, 62.4, 
200GeV in UrQMD model. 



analysis are also used in the centrality definition, it intro- 
duces an auto-correlation effect in the final results. So, 
to avoid the auto-correlation effect, we should exclude 
the corresponding proton/ anti-protons from the central- 
ity definition. 



VII. STATISTICAL ERROR ESTIMATION 

As the moment analysis is a statistics hungry study, 
the error estimation is crucial to extract physics infor- 
mation from the limited experimental data. We will 
present several statistics methods (Delta theorem. Boot- 
strap [13] , Sub-group) of error estimation in the moment 
analysis and discuss their comparisons through a Monte 
Carlo simulation. For simplicity, we assume the parti- 
cles are independently distributed as Poissonian distri- 
butions. Then, the difference of two independent Pois- 
son distributions is distributed as " Skellam" distribution. 
The skellam probability density distribution is 



/(fc;/xi,A*2) = e^^+''H^)'=/2/|fcj(2V^) (16) 

M2 



M = - /i2,cr 



(A*l+/i2)2/=^ 



1 

fH + /^2 



Scr 



/^i + Ai2 



^2 2 



1 



(17) 
(18) 

(19) 



To do the simulation, we set the two mean values of 
the Skellam distributions as /zi = 4.11, /Z2 = 2.99, which 
are similar with the mean proton and anti-proton num- 
bers measured in central Au-|-Au collisions at ^/s,^,^^ 
200 GeV by the STAR experiment Then, we gen- 

erate random numbers as per the Skellam distribution. 
The details of Delta theorem error estimation method 
for moment analysis can be found in pTj . The bootstrap 
method (20j is based on the repeat sampling with the 
same statistics of the parent distribution and the statis- 
tical errors can be obtained as the root mean square of 
the observable from each sample. In Sub-group method, 
we randomly divide the whole sample into several sub- 
groups with same statistics and the errors are obtained 
as the root mean square of the observable from each sub- 
group. In our simulation, we set 200 bootstrap times and 
5 sub-groups. 

Fig. ll2l shows the error estimation comparison between 
Delta theorem. Bootstrap and Sub-group methods for 
KCT^ of Skellam distribution. For each method, fifty inde- 
pendent samples are sampled from the Skellam distribu- 
tion with statistics 0.01 , 0.1 and 1 million, respectively. 
The probability for the value staying within ±la around 
expectation is about 68.3% and it means error bars of 33 
out of 50 points should touch the expected value(dashed 
line at unity) in Fig.[T21 We find that the results from the 
Delta theorem and Bootstrap method show similar error 
values and satisfies the above criteria, while the random 
sub-group method over estimates the statistical errors. It 
indicates that the Delta theorem and Bootstrap error es- 
timation methods for the moment analysis are reasonable 
and can reflect the statistical properties of moments. 

According to the error formula derived from delta the- 
orem, for normal distributions, the errors for cumulants 
ratios are proportional to the standard deviation of the 
distribution as: 



error{Scr) oc 



/n 



error{ka ) oc —= 



where the n is the number of events. 



(20) 
(21) 



, where the and /i2 are the mean value of the Pois- 
son distributions, respectively, the Ik{z) is the modified 
bessel function of the first kind. Then, we can calcu- 
late various moments (M, a, S, n) and moment products 
(kct^ , Sa) products of the Skellam distribution. The re- 
lations are given below: 



VIII. SUMMARY 

We presented three measurement artifacts in the mo- 
ments analysis of net-proton distributions in heavy-ion 
collisions using the UrQMD model. These are centrality 
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FIG. 12; (Color online) Ka^ for 50 samples that independently 
and randomly generated from the Skellam distribution with 
different number of events (0.01, 0.1, 1 million). The dashed 
lines are expectations value 1. 

bin width effect, centrality resolution effect and auto- 
correlation effect. The UrQMD calculations are carried 
out for Au+Au collisions at ^/s^= 7.7, 11.5, 19.6, 27, 
39, 62.4, 200 GeV. We also discussed corresponding tech- 



niques to suppress those effects to obtain physics infor- 
mation. The centrality bin width effect and centrality 
resolution effect arc originated from the volume and its 
fluctuations. The auto-correlation effect is resulted from 
using the same particles for the moment analysis and for 
the centrality definitions. We have also discussed the er- 
ror estimation method used in moments analysis. The 
effects and techniques discussed in this paper are not 
only restricted within the net-proton analysis, but also 
applicable to other particle species, such as net-charge, 
net-strangeness. 



Acknowledgments 

The work was supported in part by the National 
Natural Science Foundation of China under grant No. 
11205067 and 11135011. CCNU-QLPL Innovation 
Fund(QLPL2011P01) and China Postdoctoral Science 
Foundation (2012M511237). BM is supported by the 
DST SwarnaJayanti project fellowship. 



References 



[1] R. V. Gavai, S. Gupta, Phys. Rev. D 78, 114503 (2008). [12 
[2] M. M. Aggarwal et al. (STAR Collaboration), arXiv: 

1007.2613. [13 
Y. Aoki et al, Nature 443 675 (2006). 

S. Ejiri, Phys. Rev. D 78 (2008) 074507 ; E. S. Bowman, [14 
J. I. Kapusta, Phys. Rev. C 79 (2009) 015202. [15^ 
Y. Aoki et al, Phys. Lett. B 643 (2006) 46; M. Cheng 
et al, Phys. Rev. D 74, 054507 (2006). [16 
[6] M. A. Stephanov, Int. J. Mod. Phys. A 20, 4387 (2005); 

B. Mohanty, Nucl. Phys. A 830, 899c (2009). [17 
Z. Fodor, S. D. Katz, Phys. Lett. B 534 (2002) 87; 
Philippe de Forcrand and O. Philipsen, Nucl. Phys. B [18 
642 (2002) 290; R. V. Gavai, S. Gupta, Phys. Rev. D 
68 (2003) 034506. [19 
[8] M. A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009). 
[9] M. Cheng et al, Phys. Rev. D 79, 074505 (2009). A. [20 
Bazavov, et al, arXiv: 1203.0784. 

S. Borsanyi et al, JHEP 1201, 138 (2011). [21 
F. Karsch and K. Redlich, Phys. Lett. B 695, 136 (2011). 



[3] 
[4] 

[5] 



[7] 



[10] 
[11] 



S. Gupta, X. Luo, B. Mohanty, H. G. Ritter, N. Xu, 
Science 332, 1525 (2011). 

Y. Hatta and M. A. Stephanov, Phys. Rev. Lett. 91, 
102003 (2003). 

X. Luo (for the STAR Collaboration), arXiv: 1210.5573. 
M. Bleicher et al, J. Phys. G: Nucl. Part. Phys. 25, 1859 
(1999). 

M. M. Aggarwal et al (STAR Collaboration), Phys. Rev. 
Lett. 105, 022302 (2010). 

X. Luo (for the STAR Collaboration), J. Phys.: Conf. 
Ser. 316, 012003 (2011). 

X. Luo (for the STAR Collaboration), Acta Phys. Pol. B 
Proc. Supp.5, 497 (2012). [arXiv: 111 1.5671]. 
Michael L. Miller al, Ann. Rev. Nucl. Part. Sci. 57, 
205 (2007). 

B. Efron,R. Tibshirani, An introduction to the bootstrap, 
Chapman & Hall (1993). 

X. Luo, J. Phys. G: Nucl. Part. Phys. 39, 025008 (2012). 



